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Comparison of Fourier and model-based estimators in 
single mode multiaxial interferometry 
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ABSTRACT 

There are several solutions to code the signal arising from optical long baseline multi- 
aperture interferometers. In this paper, we focus on the non homothetic spatial 
coding scheme (multiaxial) with the fringe pattern coded along one dimension on 
one detector (all-in-one). After describing the physical principles governing single mode 
interferometers using that sort of recombination scheme, we analyze two different ex- 
isting methods that measure the source visibility. The first technique, so-called Fourier 
estimator, consists in integrating the high frequency peak of the power spectral den- 
sity of the interferogram. The second method, so-called model-based estimator, has 
been specifically developed for the AMBER instrument of the VLTI and deals with 
directly modelling the interferogram recorded on the detector. Performances of both 
estimators are computed in terms of Signal to Noise Ratio (SNR) of the visibility, 
assuming that the interferograms are perturbed by photon and detector noises. The- 
oretical expressions of the visibility SNR are provided, validated through numerical 
computations and then compared. We show that the model-based estimator offers up 
to 5 times better performances than the Fourier one. 
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1 INTRODUCTION 

The next challenge of long baseline optical interferome- 
try is to commonly perform direct imaging of the ob- 
served source, the analogous way it is done in radio- 
interferom etry (|Hogbom| |l97j) or in infrared aperture 
masking jTuthill et all I2OO0I) . A fter the first promising 
results obtain e d with COAST dBaldwin et all 11990 : 
lYoungetaD l2000h, NPOI jHummell Il998h and 
IOTA (IMonnier et all l2004ah . such technique should 
soon move one ste p forward with the operating of the AM- 
BER instrument iPetrov et al.ll200oT) . the three beam re- 
combiner of the VLTI. From the beginning of 2005, AMBER 
will indeed take full benefit of the unique combination of the 
great sensitivity of large aperture telescopes and the spatial 
frequency coverage provided by the VLTI, even though it 
will require multiple nights of observing to be able to re- 
store consistent images (iThiebaut et al J 12001 iTatulli et all 
l2004al) . Then, in less than a decade, huge improvements are 
contemplated to be accomplished with second generation 
instruments of the VLTI that will enable snapshot imag- 
ing by using 4, 6 or ev en 8 telescopes simultaneously (e.g. 
Malb et et al I (l2004h h 
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One critical point in the design of future interferometric 
imaging instruments is the choice of the beam recombination 
scheme, which can become particularly complex, especially 
when dealing with multi apertures (N te i ^ 3) interferom- 
eters. Following the solution that has been chosen for the 
AMBER instrument, we investigate the properties of single 
mode non-homothetic spatial coding scheme (from now on 
" multiaxial" ) with all the fringes pattern in the same spatial 
dimension on the same detector (from now on "all-in-one"). 
In other words, interferograms are obtained by mixing all to- 
gether the input beams arising from the different telescopes, 
thanks to output pupils arranged along one single dimension 
(see Fig.0. We analyze the ways to estimate the source vis- 
ibility from such interferograms. 

Indeed, single mode multiaxial all-in-one recombination 
appears particularly well suited in the framework of interfer- 
ometric imaging. First it is the simplest and most compact 
way to recover informations arising from all the baselines. 
Moreover it provides a better transmission than Michelson 
recombination schemes (i.e temporal coding) since it makes 
use of less mirrors and beam splitter for the same given 
number of input pupils. And the number of pixels required 
to code the signal is also smaller, which d rives to higher 
limiting magnitudes llLebouauin et af] l20041. Furthermore, 
the remarkable spatial filtering properties of single mode 
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waveguides allow to change the phase corrugations of the 
incoming turbulent wavefront into intensity fluctuations at 
the output of the fibers. In other words, only one fraction 
of the source flux, so-called coupling coefficient which de- 
pends on the Strehl ratio of the pupil apodized by the 
fiber single mode iCoude du Foresto et al.l2000fl . remains 
in the interferogram. But the very advantageous counterpart 
is that the shape of the interferogram is entirely determinis- 
tic, that is the form of the peaks in the Fourier plane 
is fixed, and the fringe pattern fully determined by 
two free parameters, its amplitude and its phase. 

Before the advent of single mode interferometers and 
in order to ov e rcome the problem of the turbulence, 
iRoddier fc Lena! lll984l) proposed to estimate the visibility 
in the Fourier plane from the integration of the high fre- 
quency peak of the long exposure power spectral densit y of 
the interferogram. Then lCoude Du Foresto et al.l B1997I) . in 
a natural way, used the same estimator to compute the visi- 
bility arising from the FLUOR experiment, the first interfer- 
ometer making use of single mode waveguides. In multiaxial 
coding, we can furthermore take advantage of the determin- 
istic nature of the interferogram shape to perform model 
fitting techniques. It was though only recently, namely for 
the AMBER instrument, that this property was used to esti- 
mate the visibility, by fitting the fringe pattern — its phase 
and its amplitude — in the detector plane jMillour et alJ 
I2004D . 

In this paper, we recall in Section |2] the general for- 
malism of multiaxial all-in-one recombination, that is the 
equation governing the interferogram, as well as the two 
techniques currently used to estimate the visibility. As men- 
tioned above, the first technique, so-called Fourier estima- 
tor, integrates the high frequency peak of the power spec- 
tral density of the interferogram whereas the second one, 
so-called model-based estimator, directly models the inter- 
ferogram in the detector plane. For both estimators, the- 
oretical expressions of Signal to Noise Ratio (SNR) of the 
visibility are provided. In Section [31 those expressions are 
validated thanks to numerical simulations of noisy interfer- 
ograms. Then in Section 0] both estimators are compared, 
from a formal point of view and in terms of relative per- 
formances. The influence of instrumental parameters is in- 
vestigated as well, with a special emphasis regarding the 
choice of the width of the reading window of the detector. 
The presence of atmospheric piston which blurs the 
fringes has not been taken into account in this anal- 
ysis. Indeed, regardless the chosen estimator, its ef- 
fect results in an attenuation of the squared visibility 
dColavitalll999^ . and the sensitivity of both estima- 
tors to this point is the same. 

This paper is the first part of our study on 
mutiaxial all-in-one recom bination. In a second paper 
jLebouauin fc TatulliEooil) . we analyze in which way such 
recombination scheme allows to optimize the visibility SNR 
thanks to specific geometric configuration of the output 
pupils. 



2 GENERAL FORMALISM 

Figure sketches the principle of multiaxial recombination 
in waveguided interferometers. The light arising from the 
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Figure 1. Sketch of a multiaxial all-in-one single mode interfer- 
ometer. t l denotes the total "static" transmission from the i th 
telescope aperture to the optical plane A (i.e. mirrors, delay line, 
transmission of the i th fiber, ...), whereas p 1 takes into account 
the "dynamical" transmission, that is the coupling coefficient of 
the i th fiber. i(a) is the interferogram. p x (a) and pi (a) are the 
photometric channels. The photometric fluxes F l , F J and the 
coherent flux F^ 3 are defined in the optical plane A. a 1 (a) and 
6* (a) are respectively the detected beams in the interferometric 
and photometric channels. In other words, a' (a) and b l (a) are 
the transmission factor between the photometric flux and the in- 
terferometric and photometric channels, respectively. Note that 
the definition of a 1 (a) and b l (a) includes the transmission of the 
beam splitter. 

Intensity 

1 -5" Npattem 



Nfri 




Figure 2. Image of a typical interferogram at the output of mul- 
tiaxial all-in-one recombiner. Here is shown a 2-beam recombina- 
tion with instrumental contrast set to 1. The fringes are weighted 
by the intensity mode pattern (dashed line). Npattem in- 
dicates the width of the detector reading window, the 
width of one mode pattern being defined as the interval 
between the top of the Gaussian mode and where the 
amplitude of the Gaussian mode has decreased by a fac- 
tor exp(— 7r). This definition has been chosen so that this 
width corresponds to half of the first lobe of the Airy 
pattern in case of pupils which are not weighted by the 
first mode of single-mode fibers. Nf ri refers to the num- 
ber of fringes in one mode pattern. This term is fixed by the 
distance between the output pupils. 

i th telescope is filtered by a single mode fiber to convert 
phase fluctuations of the corrugated wavefront into inten- 
sity fluctuations. The fraction of light p 1 entering the fiber 
is called the coupling coefficient jShaklan fc R,oddiedll98sll 
and depends on the strehl ra tio of the apodized pupil 
llCoude du Foresto et al.|l2000t) . Making use of a beam split- 
ter, one part of the light is selected to estimate the photome- 
try, thanks to dedicated photometric channels. The remain- 
ing part of the light is recombined with the beam coming 
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for the j telescope to form fringes. The coding frequency 
of the fringes f 13 is fixed by the separation of the output 
pupils, which are arranged along one dimension. 

When the only i th beam is lighted, the signal recorded 
on the interferometric channel is the photometric flux F l 
spread on the intensity mode pattern a 1 (a), that is the 
diffraction pattern of the i th output pupil weighted by the 
single mode of the fiber, a is the angular variable in the im- 
age plane. F* results in the source photon flux N attenuated 
by the total transmission of the instrument, i.e. the product 
of the "static" transmission t z and the coupling coefficient 
p 1 of the single mode fiber: 



F l = Nt l p l 



(1) 



When beams i and j are lighted simultaneously, the 
coherent addition of both beams results in an interferomet- 
ric component superimposed to the photometric continuum. 
The interferometric part, that is the fringes, arises from the 
amplitude modulation of the coherent flux Fl 3 at the coding 
frequency f tJ . The coherent flux is the geometrical product 
of the photometric fluxes, weighted by the visibility: 



Fl 3 = 2NV¥Uy / p i pjV i - 



(2) 



wher e V' 3 e^ 3 is the complex modal visibility llMeee et alJ 



2001) and ^ takes into account a potential differential at- 
mospheric piston. Note that strictly speaking the modal vis- 
ibility is not the source visibility. Rigorously, the modal 
visibility depends on the convolution between the 
source visibility and the telescope transfer function 
which is atmosphere dependent in the optical range. 
As such the modal visibility is biased both by the ge- 
ometric antenna-lobe effect (the object is multiplied 
by the telescope point spread function, as commonly 
known in radio-astronomy), and by the turbulence. 
When the object is unresolved by one single tele- 
scope however, the modal visibili ty can be f airly 
approximated by the object one (Tat ulli fc C helli 
2005) and its estimation is robust (stable to the level 



of 1% or less) to a change of atmospheric during 
the calibration process (Tatul li et al.ll2004bf) . In any 
case, a further study of the relationship between the 
modal visibility and the source visibility is beyond 
the scope of this paper, and further informations can 
be found in papers mentioned above. Here we consider 
our observable to be the modal visibility. 

Such an analysis can be done for each pair of telescopes 
available in the interferometer. As a result, the interferogram 
recorded on the detector can be written in the general form: 



z(a) = J»)F l + 

i 

NteZ 

£ Va«(a)oi(o)Cg(a)Re [i^V (2 ™ rJ+ ^ (a)+ *^} 

where <f>l 3 (a) is the instrumental phase taking into account 
possible misalignment and/or differential phase between the 
beams a z (a) and a 3 (a). Cg and <& l g(a) are respectively 
the loss of contrast and the phase shift due to dif- 
fusion and polarization effects 1 , which may not be 

1 assuming a non-polarized incoming light 



homogeneous along the fringe pattern. Fig. Ogives an 
example of a multiaxial all-in-one interferogram in the two 
telescope case. 

Thanks to the photometric channels, the number of 
photoevents p l (a) coming from each telescope can be es- 
timated independently: 

p l (a) =F l b\a) (4) 

where b l (a) is the detected beam in the i th photometric 
channel. 

We can notice from Eq.'s (0 and |J5) that the estimator 

of the squared modal visibility |V lJ | 2 results in the ratio be- 
tween the squared coherent flux and the photometric fluxes. 
Using the previous definitions, we can set a generic form of 
the estimator as following: 

< IK 3 1 2 



\V i3 \ 2 



> 



(•») 



< 4F*FJ > 

Note that < \F* 3 \ 2 > is computed instead of < Fl 3 > be- 
cause in absence of fringe tracking the random atmospheric 
differential piston (j>p totally blur the coherent signal. It now 
remains to estimate F 1 and Fl 3 from the interferogram. 

2.1 Fourier estimator: integrating the power 
spectral density 

In the Fourier space, the interferogram defined by Eq. JHJ 
takes the form of the sum of photometric and interferomet- 
ric peaks. The photometric peaks are centered at the zero 
spatial frequency whereas the interferometric peaks MV(f) 
are located at their respective spatial coding fre- 
quency their counterpart Mj!(f) being in the 
negative spatial frequency domain. This method in 
the Fourier space requires that (i) the photometric 
peaks and the interferometric peaks are not over- 
lapping, and, (ii) the high frequency peaks are not 
overlapping between each other. If these conditions 
are fulfilled, the squared coherent flux can be estimated 
by comput ing the integral of the power spect ral density 
\M^(f)\ 2 dRoddier fc Lenalll984l : IConarJ Il994l) . over its 
frequency support, that is [f %3 — D/X, f lJ +D/X] where 
D is the diameter of the output pupil. From the defi- 
nition of the coherent flux and using the Parseval equality, 
it comes: 

The photometric flux is easily computed from the photomet- 
ric channel (see Eq. Q): 



F l \ 6 i (a)da 



Then the estimation of the fringe contrast C 13 writes: 



with 



< P*- Pi > 



n 2_jC l i {a)a\a)a 3 \a)da 



/ ¥(a)da J b 3 (a)da 



(7) 



(8) 



(9) 



being the instrumental contrast of the recombiner that de- 
pends on the contrast loss due to polarization effects, 
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on the alignment of the beams a 1 (a) and a 3 (a), and on the 
flux ratio between the interferometric and the photometric 
channels. Note that the power spectral density of the in- 
terferogram ha s to be prop erly unbiased from photon and 



terrerogram nas to be propc 
detector noise dPerrinll200al . 



= [CI 



R 
I 

F 



where the matrix C takes the detailed form: 



(14) 



2.2 Model-based estimator: modelling the 
interferogram 

The model-based estimator has been introduced for the first 
time i n the data reductio n process of the AMBER instru- 
ment jMillour et al.l2 0041. It consists in modelling the inter- 
ferogram thanks to a priori knowledges of the instrument. 
The purpose of such a signal processing is twofold: (i) to 
develop optimized algorithms in terms of performances of 
the instrument, that is the SNR of the visibility; and (ii) 
on the contrary of the Fourier estimator, to authorize high 
frequency peak overlapping when dealing with multi-beam 
[Ntei > 3) recombination, thus allowing to code the interfer- 
ogram on less pixels. This second point is beyond the scope 
of this paper. Let just mention here that plainly choosing 
the different coding frequencies is crucial to design opti- 
mized multi-beam recombiner making use of integrated op- 
tics. This is especially tr ue in the case of imagi ng instru- 
ments such as VITRUV l|Lebouauin et alJ 120041) that are 
recombining 4 beams or more, as it is shown in our second 
paper on the subject jLebouauin fc TatulliEooi) . 

A full descripti on of this estimator can be found in 
iMillour et alJ l|2004h . We only recall here the basics prin- 
ciples. To model the signal on the detector, Eq. has to 
be rewritten in its sampled version, where k stands for the 
pixel number, between 1 and N p i x : 



i i<j 



I! + 4 i ' i) J r " 



(10) 



with 



and 



d[ t ' 3) and F 3 being the quadratic counterpart of c k and 
R' J respectively. c^' J ' and dl are called the carrying waves 
of the signal at the coding frequency f lJ , since they " carry" 
(in terms of amplitude modulation) R 13 and F J , which are 
directly linked to the complex coherent flux. Furthermore, 
the photometric fluxes are still computed from the photo- 
metric channels (see Eq. 0): 



C%[k) 



\lY. k cf{k)aia 



3(27ra fc / y +cbl 3 (k) + $ i i(k)) 
(11) 



(12) 



(M) 



(13) 



F % and Fl 3 are then jointly estimated from the photometry 
[p 1 ) and the interferogram (i k ) by resolving a set of [Npix + 
N te i) linear equations with (2Nt + N te i) unknowns [Nb being 
the number of pairs of telescopes, i.e. Nb = Nt e i[N te i — l)/2): 



( ■■ 







Aid) 



a i 











\ .. .. .. £6^< / 

(15) 

The matrix C entirely characterizes the instrument. It de- 
pends on the shape of the detected beams a\, b\ and on the 
carrying waves c k ' 3 \ d k % that hold informations about the 

interferometric beam y a 'k a k> l ne coding frequencies f 3 and 
the instrumental differential phases cf> 13 , and the polarization 
state within Cg and "1?^ . Such quantities can be calibrated 
in laboratory, hence they are assumed to be pe rfectly known. 
The ca libration procedure is fully described bv lMillour et alJ 
(2004). C has to be inverted in order to solve the system. In 
the AMBER experiment, the generalized inverse of C has 
been called the Pixel To Visibility Matrix (P2VM), since it 
enables to compute the visibility of the fringes from the mea- 
surements on the detector. The estimation of the contrast 
writes: 



\cw 



< R 13 ' 1 > + < I ij2 > 



\vv\ 2 .c; (16) 



Cf, is still the squared instrumental contrast: 



C r 



Z h Cri\k)aiai 



(17) 



with the same definition than in Eq. ©. Note that the quan- 

■ '2 '2 

tity R 13 + F 3 has to be properly unbiased, like the power 
spectral density in the Fourier plane. 



2.3 SNR of the modal visibility 

Using second order expansion of lPapoulisI (|W84) , we derive 
from Eq. the relative error (i.e. the inverse of the SNR) 
of the squared visibility: 



£ 2 (\V 



a 2 [\F^\ 2 ) 



(F 3 



+ 2 



Cov(F\F 3 ) 
Fi-Fi 



F r Fi~ 

Cov(\F?\ 2 ,F l F 3 ) 



J?V 12 F i F j 



(18) 



The main difference between both approaches lies in the 
following remark: in the case of the Fourier estimator, the 
coherent and photometric fluxes are directly estimated from 
the measurements, each independently, whereas in the case 
of the model-based estimator, the coherent and photometric 
fluxes are jointly reconstructed from the measurements by 
way of computation of the P2VM matrix, and are therefore 
correlated. As a result we have the following situations: 
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Fourier estimator Model-based estimator 



Cov(|F« | 2 , F^-F 3 ') = Cov^ 3 \ 2 ,F'F 3 ) / 
Cov(F l ,F J ) = Cov(F\F J ) / 

Detailed computation of Eq. (|18^ are given in Appendix 
lAl for both estimators, assuming that the interferogram is 
corrupted by photon and detector noise. Atm ospheric noise 
is negle cted here since it has been shown in iTatulli et alJ 
(2004b) that in presence of modal filtering — and on the 
contrary of multispeckle interferometry l|Goodmanl 
ll985lT) — , speckle noise is the dominant noise in the 
case of very bright sources (negative magnitudes) 
only, and is therefore marginally relevant. Further- 
more, we did not take into account in these com- 
putations the effect of atmospheric piston since it 
only results in the attenuation of the squared visibil- 
ity, the sensitivity of both estimators to this specific 
point being the same. 

In next section, we propose to simulate the estimation 
of the modal visibility from multiaxial recombination and to 
validate our theoretical calculations. Then we compare the 
performances of both estimators. 



3 VALIDATION OF THE THEORETICAL 
EXPRESSIONS 

In order to validate our theoretical expressions derived in the 
previous section, we perform statistical simulations of noisy 
interferograms. For sake of simplicity we assume from now 
on a two telescope interferometer. Note however that the 
validity of our theoretical approach has been also checked for 
increasing number of telescopes. Moreover, a deep analysis 
of multi-beam (55 3) r ecombination is proposed in our second 
paper on the subject jLebouauin fc TatullilEooil . 

Following the formalism of Section D2J the interferogram 
arising from multiaxial all-in-one recombination is entirely 
defined by the following parameters: 

(i) The shape of a(a), which is assumed to be the same 
for each beam. It arises from the inverse Fourier Transform 
of the auto-correlation of the output pupil of diameter D. To 
take into account the weighting of the single mode 
fiber we assumed a Gaussian shape with a pupil 
stop of width D (that is a convolution by a Bessel 
function in the detector plane). As a result, the low 
and high frequency peaks have truncated Gaussian 
shapes with a base width of 2D/X. We consider that the 
photometry is recorded on one pixel. 

(ii) The number of fringes Nf r i in the interference pattern 
for the lowest coding frequency. It is defined as the distance 
between the closest output pupils in D/X units (see Fig.[5Jl. 
We choose here Nj r i — 2. 

(iii) The number of pixels per fringe N p f to code the inter- 
ferogram. It must be chosen such that it fulfills the Shannon 
criteria for the highest frequency coding f ma x of the carrying 
waves. If it is written under the form f m ax = /3D/X, then 
the number of pixels must verify: N p f ^ 2/3. For the two 
telescope case considered here, we arbitrarily set N p f — 2.5. 

(iv) The width of the detector reading window N pa ttem 
which fixes the total number of fringes taken into account in 
the interferogram and hence the total number of pixels N p i x 



read on the detector. We impose in this section N pat tern = 2. 
This choice means that the detector reading window is two 
mode patterns wide, exactly as if we only would con- 
sider the fringes in the first lobe of the diffraction 
pattern, in the case of an Airy disk. Such choice seems 
reasonable at first thought since outside this lobe the inter- 
ferogram is severely attenuated (as one can notice in Fig. [5J . 
Nevertheless, a deeper analysis of this specific point shows 
that such a parameter is a key issue, as it will be discussed 
in Section l4~2l 

(v) The fraction of flux going into the photometric chan- 
nels (that is the transmission of the beam splitter). We as- 
sume here that the beam-splitter selects 30% of the flux for 
the photometry. 

For a given source magnitude, the number of photoevents 
occurring on each pixel of the interferometric and the pho- 
tometric channels are computed following Eq.'s and @ 
respectively, assuming photon noise and additive detector 
noise with a = 15e~/pix. Such a procedure is then repeated 
until we obtain a sample of 1000 data sets, which is large 
enough to perform statistics. For both estimators, we com- 
pute theoretical and statistical mean value and Signal to 

Noise Ratio of the modal visibility |V ,J '| 2 , thanks to Eq.'s 
JIjJ and 11811 . Fig. shows the results of our computations, 
for both methods (formal and simulated) and for both es- 
timators (P2VM and Fourier), setting the true value of the 
modal visibility to 0.5. Theoretical calculations and numer- 
ical simulations are in excellent agreement, both for the es- 
timated visibility and the SNR. This study validates the 
theoretical expressions of both estimators as well as their 
respective theoretical SNR. 



4 DISCUSSION 

4.1 Estimator relative performances 

Although both estimators arise from the same formal defini- 
tion of Eq. JKJ , they exhibits fundamental conceptual differ- 
ences. Obviously, both techniques presents the same instru- 
mental contrast (see Eq.'s @ and (1171 respectively), which 
is not surprising since the instrumental design is strictly the 
same, as well as the same estimation of the photometry (see 
Eq.'s (|7J and 11311 respectively). However the very difference 
lies in the computation of the coherent flux. 

First, in the Fourier case, the coherent flux results in 
a second order (i.e. quadratic) estimation, that is | i 7 ^ J 1 2 
is directly estimated. On the contrary, the model-based 
computation is equivalent to a first order estimation, that 
is the complex quantity Fl 1 is calculated, and then the 
squared modulus is taken to get rid of the atmospheric dif- 
ferential piston. This latter method is equivalent to 
fit the complex Fourier Transform of the interfer- 
ogram. A second order estimator based on mod- 
elling would have been the fitting of the power spec- 
tral density itself, or equivalently, the fitting of the 
auto-correlation of the interferogram in the detec- 
tor plane. Performing first order estimation is above 
all interesting because it allows to separate infor- 
mations (i.e. the visibilities for each baseline) be- 
fore computing the modulus square. As a result, 
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Figure 3. Theoretical (solid lines) and simulated (cross) values of the modal visibility mean value (top) and SNR (bottom), for the 
P2VM (left) and the Fourier (right) estimators. Plots are shown as a function of the number of photoevents per pixel in the interferogram. 
Detector noise has been set to a = 15e — /pix. The vertical dashed line shows the limit between the detector and the photon noise regime. 
For the estimation of the visibility arising from numerical simulations, we plot the statistical error bars due to limited number of samples 
(1000 data sets), that is the dispersion of the 1000 estimated modal visibility, divided by the square root of the number of samples. 



this method prevents "cross-talk" between the base- 
lines even if the peaks are partially superimposed, 
which is particularly worthwhile in the case of vis- 
ibility estimation from multi-beam (Js 3) interfer- 
ometers, as develope d in our secon d paper on the 
subject (Lebouauin &: Tatulli 2005). We also infer 
that a first oder estimation drive to better — or at 
worst identical — performances than quadratic one, 
though a thorough analysis of this point, which is 
beyond the scope of this paper, remains to be done. 
Moreover and above all, the model-based algorithm, thanks 
to the P2VM calibration matrix, takes entire benefit of the 
knowledge of the instrument whereas the Fourier estimator 
does not 1 . When making use of the model-based estimator, 
the shape of the interferogram is perfectly known, precisely 

its envelope (the diffraction pattern <J a^aj.) as well as its 
coding frequency and its instrumental phase. These a priori 
informations are gathered in the matrix C defined in Sec- 
tion |2i3 And clearly, introducing perfectly known (i.e true 
and un-noisy) a priori in the data reduction procedures can 



1 although some approximate assumptions about the peak posi- 
tion have to be made to set the integration interval. 



only improve the performances of the corresponding estima- 
tor. These two remarks, and particularly the second point, 
explain why the model-based algorithm leads to better per- 
formances than the Fourier one, as it is illustrated in Figure 
0] Note that the SNR improvement is all the more important 
than the visibility is high. In the two telescopes case, a gain 
of a factor 3 to 5 can be achieved at best, for unresolved 
sources. 

However this analysis assumes perfect calibration of 
the instrument. It means that the calibration matrix must 
be both perfectly stable in time and very precise, that is 
recorded with a SNR much higher than the SNR of the in- 
terferograms. If the instrument is not stable between the 
calibration procedures and the observations, the P2VM will 
drift and as a result, the estimated visibilities will be biased. 
And if the calibration is not precise enough, it will be the 
limiting factor of the visibility SNR. In the case of the AM- 
BER instrument, the calibration procedure is quite complex 
and it can require a consequent integration time (several 
minutes) to get a useful ad precise calibration. More gener- 
ally, the time and the way to calibrate an instrument severely 
depends on its stability and on its complexity. Ambitious 
designs su ch as the "silicon v- gro ove array" of the MIRC re- 
combiner jMonnier et al . 2004b), or recombination schemes 
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nb. of coherent y / pix 

Figure 4. Ratio between the SNR of the Fourier estimator and 
the SNR of the model-based one, as a function of the number of 
photoevents per pixel in the interferogram. Detector noise is still 
a = 15e _ /pix. The curves are plotted for 3 types of sources: fully 
resolved (V = 0.05, solid line), moderately resolved (V = 0.5, 
dashed line) and unresolved (V = 0.95, dotted line). 

making use of integra ted optics chips feerger et al]l2003t 
iLebouauin et al.l2004) should drive to drastic improvements 
on this specific point. 

4.2 Influence of instrumental parameters 

Previous analysis has been done with a given configura- 
tion of the instrument. One last point to investigate is how 
both estimator behave, one compare to the other, when the 
parameters governing the interferogram are varying. Obvi- 
ously, modifying the number of pixel per fringe or even the 
detector noise level will result in similar changes for both 
estimators, that is the slope of the SNR in the detector 
noise regime. As well, changing the coding frequency, which 
only defines the position of the interferometric peak in the 
Fourier space, will lead to equivalent modifications of both 
estimators performances, at least as far as the high frequency 
peaks are separable 2 . At last, choosing the optimized area on 
which the interferogram gives valuable informations without 
adding too much detector noise is a crucial point. Since the 
model-based algorithm takes into account the shape of the 
interferogram whereas the Fourier estimator does not, the 
compromise to find is not the same in both cases. It means 
the response to a change of the detector reading window 
will differ with regards to the chosen estimator. The effects 
of this parameter are investigated here. 

In a multiaxial combination, one has the choice of the 
limits of the reading window on the detector, i.e of the num- 
ber of pixels to consider. And the largest the window, the 
more signal you integrate, but the more detector noise you 
record too. In Figure |3J we show the evolution of the SNR 
ratio of the modal visibility as a function of the width of 
the window on the detector (here defined in fraction of in- 
tensity mode pattern). The entire first lobe of the fringe 

2 when N tt ,t ^ 3, see lLebouauin fc Tatulhl <2005l) 



pattern contains Nf r i — 8 fringes, with N p f — 4 pixels per 
fringe. All the other parameters of the instrument are kept 
unchanged. 

Fourier estimator: In the photon poor regime, at the 
detection limit of the instrument, the SNR shows a maxi- 
mum when the width of the reading window is about half 
of the first lobe of the envelope. Beyond that point, the pix- 
els have very small individual SNR and only bring a noise 
contribution in the estimation of the visibility. In the pho- 
ton rich regime, the SNR reaches its maximum for the same 
width of the reading window but exhibits a plateau as the 
width increases. This behavior stands as long as the SNR 
of each individual pixel is dominated by the photon noise. 
Then the SNR of the visibility starts to decrease. The width 
of the plateau depends on the incoming flux of the source 
and is all the more large than the source is bright. 

Model-based estimator: As for the Fourier estimator, the 
shape of the SNR is linked to the shape of the fringe pat- 
tern. But in this case, the SNR is increasing with the width 
of the detector window and does not exhibit a maximum (in 
other word, the optimal is found for an "infinite" width). 
As a matter of fact, thanks to the generalized inverse of the 
matrix C that takes into account the shape of the interfer- 
ogram, each pixel contribution is weighted by its individual 
SNR. So the pixels with bad SNR (due to the envelope or 
fringe modulation) are "removed" from the reconstruction 
and do not introduce noise in the estimation of the visibility. 
Nevertheless, we can see that the slope of the SNR becomes 
almost flat from a detector width of about one intensity 
mode pattern. 

This analysis shows that, in the framework of inter- 
ferometric observations making use of multiaxial all-in-one 
recombination, and in the case of bright sources, it is worth 
the effort to integrate the interferogram on the entire 
lobe in order to optimize the SNR of the visibility. This 
statement stands for both estimators. When observing faint 
sources, that is when reaching the limiting magnitude of the 
instrument, and in the specific case of the Fourier estima- 
tor, performances are slightly improved when reducing the 
width of the reading window to half of the lobe, although 
the gain on the SNR never exceed a factor of 2. 



5 CONCLUSION 

In this paper, we have developed the theoretical formalism 
that allows to model single mode interferometers using mul- 
tiaxial all-in-one-coding, from the signal processing point of 
view. From this formalism, two estimators of the visibility 
have been analyzed. The first one consist in using the clas- 
sical integration of the power spectral density of the inter- 
ferogram in the Fourier plane whereas the second one deals 
with modeling the interferogram in the detector plane, as it 
has been chosen for the AMBER experiment. Performances 
of such estimators have been computed. Considering photon 
and detector noises, theoretical expression of the SNR of the 
visibility have been recalled for the Fourier estimator, and 
derived for the first time in the case of model-based estima- 
tors. These expressions have been validated through numeri- 
cal simulations and then compared. We have shown that the 
second technique offers optimal performances since it makes 
full use of knowledges about the instrument, especially the 
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Figure 5. SNR of the modal visibility as a function of the detector reading window, in fraction of intensity mode pattern, for photon 
poor (left, here 17/pix ) and photon rich (right, here 10 6 7/pix) regimes. Results are displayed for marginally resolved (V = 0.9, top) 
and fully resolved (V = 0.1, bottom) sources. The instrumental parameters are: Nf ri = 4 fringes per beam, N p f = 4 pixels per fringe, 
and a detector noise of a = 15e — /pix. The vertical dash lines correspond to recording one entire mode pattern (see Fig. ISt. 



shape of the interferogram. In the two telescopes case which 
has been emphasized in this paper, we have demonstrated 
that the model-based estimator enables at best to come over 
a factor of 5 of the visibility SNR, compare to the Fourier 
one. Finally we have addressed the question of the width 
of the reading window of the detector. This point is indeed 
a crucial issue when dealing with multiaxial recombination. 
We have shown that, regardless of the chosen estimator, in- 
tegrating the entire lobe of the intensity mode pattern offers 
optimized performances. 
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APPENDIX A: THEORETICAL SNR OF THE 
VISIBILITY 

We assume that the interferogram as well as the photometric 
outputs are corrupted by photon (Poisson) noise and addi- 
tive Gaussian noise of variance a 2 . 



Al Generic expression 

The estimator of the square visibility can be express in a 
generic form: 



\V^\ 2 oc 



< \F^\ 2 > 

< F % Fi > 



(Al) 



where Fl 3 the coherent flux at the frequency f v , an F\ F' J 
are the photometric fluxes. 

The relative error S(\V ij \ 2 ) on the square visibility is 
then given bv lPapoulisI (Gi?84): 



£ 2 (\V 



a 2 (\F^\ 2 ) , <r 2 (F l ) , a 2 (F 3 ) 



+ 



+ 



F' 



; n Cov(F',^) 2 Covgg^) 



F' Fi 



Theoretical expression of each term of previous equation is 
now given for both estimators. 



Hence the expected value and the error on the coherent flux 
writes: 



|F c lJ | 2 



^(IF^I 2 ) oc 



k 

^> 2 (!AT(/ fc )| 2 ) + 



(A4) 



^^CovdM 1 ^/,)! 2 ,^'^)! 2 ) (A5) 

fe l^k 



The statistics of the spectral de nsity of an interf erogram 
have already been computed by |Goodm^J^j85[) in the 
case of photon noise and completed bv lTatulli et ajj (2004b) 
with detector and atmospheric noise. We recall the results 
here, without taking into account the atmospheric noise (i.e. 
speckle noise): 



|M*i(/ fc )| 2 = N \i(f k )\ 2 + N + N plx a 2 det (A6) 

where we recognize the bias part due to photon noise 
(JV, [Goodman (1985)) and additive Gaussian noise 
(A^gj^. lTatulli et alj d2004bh ). 

a 2 (\M l] (f k )\ 2 ) = 2iV 3 |?(/ fc )| 2 +4iV 2 |i(/ fc )| 2 + iV 2 + 



TVjLa 4 + 3N plx a 4 + 2N plx a 2 N + 

2TvF 2 f 



2N pix a z N \i(f k )\ 



(A7) 



Cov(|M lJ (/ fc )| 2 ,M lJ (/ ; )| 2 ) = 27V 3 Re \i(f k )l*(fi)i*(h - ft 
2N 3 Re [<(/i,)i(/,)i*(/* + fl)] + 2iV 2 R/ fe )| 2 + 27V 2 |T(/ i )| 2 
N 2 \i(f k - fi)\ 2 + N 2 \i{fk + fi)\ 2 + N + SN ptx a 4 

where i(f) is the normalized spectral density (such as i(0) = 
1), that is: 



Kf) 



V(f)g(f) 



(A9) 
(A10) 



for the photometric and the interferometric peaks 
respectively, N te i being the number of telescopes, 
and g(f) being the normalized Fourier Transform of 
the beam a(a) (assuming, the shame shape for the 
whole beams). Furthermore, we have 



a 2 {F l ) =F> + aj e 



(All) 



Finally, since the coherent flux and each photometric flux 
are estimated independently, we have: 



Cov(|F ( 



ij\2 F i F j 







Cov(F\F J ) = 



(A12) 



(A13) 



A2 SNR for the Fourier estimator 



A3 SNR for the P2VM estimator 



The coherent flux is linked to the spectral density of the in- 
terferogram by the following relationship (see Eq. (|SJ , writ- 
ten in its sampled form): 



|F^| 2 cx^|M^(/ fc )| 2 



(A3) 



We recall that the real and imaginary part of the weighted 
complex visibility are defined by the system of equations: 



i 

P 



= [C]. 



R 
I 
F 



(A14) 
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if we call M = mu, k S [l--N P i X +N te i] the vector resulting in 
the concatenation of the interferogram i and the photometry 
P, we can write: 

N ■ 

R l] = £ a j m k (A15) 
k 

I 13 = £ (A16) 



F* = E ftm* (A17) 
fc 

where ^ and /JJ are the coefficients of the P2VM matrix. 
Hence it comes: 

| F ii|2 = + jiP =Y J YM^l+^Q i ]m k m l (A18) 

F'f 3 ^^^m t m, (A19) 
fe i 

Here, the covariance between the coherent flux and the pho- 
tometric fluxes, as well as the covariance between the pho- 
tometric fluxes have to be taken into account. For sake of 
simplicity, Eq. I|A2(I can be rewritten: 

£ a (|V«| a ) = + ° 2 ( F ' FJ ) _ 2 Cov(\F?\ 2 ,F*F>) 

ipij |2 2 F i Fi 2 \F^\ 2 F i Fi 

(A20) 

It now remains to compute all the terms knowing that: 

a 2 {\F l Jf) = JfJ^-JfJ 1 ^ 2 (A21) 

a 2 {F i F j ) = F i2 Fi 2 -WWi 2 (A22) 

CovflF^FVF 3 ') = IF^p^FJ - \Fl 3 \ 2 WW (A23) 

To lighten the calculations we introduce the variable 7 such 
that = CkCi 3 + Ck Q 3 ■ Then we can compute the second 
order statistics of the square coherent flux and the photo- 
metric fluxes: 

\Fc | 2 = E E ik{mmi = E ik k ml + E E 7%™* m 

k l k k Ij^k 

(A24) 



Table Al. Coefficients of the fourth order statistics of the coher- 
ent and photometric fluxes. 
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computed knowing the statistics of rrifc, that is: 

ml = Wk 2 + mk + a 2 (A27) 
p 2 = p 2 + p + N pho tpixV 2 (A28) 

Nphotpix being the number of pixels to code the pho- 
tometric outputs. The third and fourth moments m| 
and m| are derived from first and second order ones 
assuming Gaussian statistics for sake of simplicity. 



W = £ piPlm 2 + E E ™i (A25) 

fe fe Ijtk 

Then the fourth order statistics F 1 = F c i3 | 4 , F i2 Fi 2 , 
\F% 3 \*F*Fi can be described by the generic equation: 

f~ 4 = E^+EE^ 7777 ^ 

* fe 1^*; 

fe i^fe i^fc^n 

EE E E "fetaa 777 ^ 7777 ^ ^ ( A26 ) 

, (1) (2,1,1) (2,2) (3) j (4) ■ ■ , 

Where "feino- "fein > Q fei . a kl 7 and «* are S 1Ven 111 t& ~ 

ble lAll in each specific case. Previous equations can be 



